The phase structure of a chiral model with dilatons in hot and dense matter 
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We explore the phase structure of a chiral model of constituent quarks and gluons implementing 
scale symmetry breaking at finite temperature and chemical potential. In this model the chiral 
dynamics is intimately linked to the trace anomaly saturated by a dilaton field. The thermodynamics 
is governed by two condensates, thermal expectation values of sigma and dilaton fields, which are 
the order parameters responsible for the phase transitions associated with the chiral and scale 
symmetries. Within the mean field approximation, we find that increasing temperature a system 
experiences a chiral phase transition and then a first-order phase transition of partial scale symmetry 
restoration characterized by a melting gluon-condensate takes place at a higher temperature. There 
exists a region at finite chemical potential where the scale symmetry remains dynamically broken 
while the chiral symmetry is restored. We also give a brief discussion on the sigma-meson mass 
constrained from Lattice QCD. 

PACS numbers: 12.39.Fe,12.39.Mk,12.38.Mh 
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1. INTRODUCTION 

Effective theories of strongly interacting matter are 
expected to capture non-perturbative aspects of QCD 
in low-energy domain. They are constructed based on 
global symmetries of QCD Lagrangian and their breaking 
pattern. In the limit of massless quarks the Lagrangian 
possesses the chiral symmetry and scale invariance, both 
of which are dynamically broken in the physical vacuum 
due to the strong interaction. The QCD trace anomaly 
signals the emergence of a scale at the quantum level from 
the theory without any dimension- full parameters [l|. 
Thus spontaneous chiral symmetry breaking, which gives 
rise to a nucleon mass, and the trace anomaly are closely 
linked to each other [2| and dynamical scales in hadronic 
systems are considered to originate from them. How they 
behave under extreme conditions such as high tempera- 
ture and density is one of the main issues in QCD [3j. 

The trace anomaly has been implemented in a chiral 
Lagrangian by introducing a dilaton (or glueball) field 
representing the gluon condensate (G^G^) Q- Ther- 
modynamics of the dilatons at finite temperature and 
density has also been explored and the deconfinement 
phase transition was studied [5|- Incorporating the QCD 
scaling properties into a non-linear chiral Lagrangian, 
the in-medium scaling associated with chiral symmetry 
restoration, BR scaling [y), was introduced and some re- 
lated works have been carried out @, H[. Besides, along 
with the Lattice QCD computations, pure gluon dynam- 
ics at finite temperature has been formulated in several 
approaches |9l-[l3j. 

In this paper we introduce a model of constituent 
quarks and gluons implementing chiral and scale invari- 
ance in such a way that the model mimics the non- 
perturbative nature of QCD in low energies. We will 
explore the thermodynamics and constrain the sigma me- 
son mass utilizing the QCD trace anomaly extracted from 
Lattice QCD pjj]. Imposing field theoretical require- 



ments on the anomaly matching, we will give a suggestive 
phase diagram of QCD. 



2. A TOY MODEL 

In this section we briefly introduce our model for con- 
stituent quarks and gluons restricting to a system with 
two flavors. 

Scale invariance is implemented in a linear sigma model 
via the following Lagrangian 



C = qiftq + Gsq (c + if ■ ff) q 



+ o -v„-v x , 
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(2.1) 



where G$ is the scalar coupling constant and B is the bag 
constant. All other notations follow the standard linear 
sigma model. We assume that the constituent gluons be- 
come massive due to the non- vanishing gluon condensate, 
(x) 0- This is achieved by introducing the Lagrangian 
for the constituent gluon field A M , 



C A =- -A^ v + -G\ [ ± ) A»A» , (2.2) 



* x There are some uncertainties on introducing \ m the explicit 
breaking term. See e.g. [^. Il4ll . This does not change our results. 
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with the field strength tensor A^ u — d^A u — d v A^ and 
the coupling constant Ga to the dilaton field. The full 
Lagrangian is thus given by 



C + jCa- 



(2.3) 



Here we assume that the quarks have no direct coupling 
to the gluons since the interaction between the quarks 
and gauge fields is embedded in Gs and Ga- 

Applying the mean field approximation, one finds the 
thermodynamic potential by performing the path inte- 
gration over the quark and gluon fields: 

n = n + n A + v a + v v + \B, 



&A = -1A 



d 3 p 

(2tt) 3 
d 3 p 



(27T) 



Tpn(l-n 9 ) + ln(l-n g )] , 
3 Tln(l + n A ), (2.4) 



with the degeneracy factors for quarks 7 9 = 2NfN c = 12 
and for gluons ja = 2(N 2 — 1) = 16. A constant term is 
added so that Sl = 0atT = /i = 0. The effective masses 
of the quasi-particles are defined by 



M q = G s a, M A = G A —. 

Xo 

The thermal distribution functions are given by 
1 1 



(2.5) 
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• e (E q -n)/T + 1 > ll 1 e (E q+l i)/T + i ' 
n * = e E A /T _ 1 ' ( 2 - 6 ) 

with the quasi-particle energies E q — ^J\p\ 2 + M 2 and 

0, leads to the 

X 



Ea = VW + M\. 

The stationary condition, |^ = |p 
following coupled gap equations: 
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The mesonic parameters A and e are related with the 
sigma and pion masses and the pion decay constant via 



A 



Vi 



(2.9) 



where the vacuum sigma expectation value is Co — fn- In 
the following calculation we will use = 138 MeV and 



/,r = 93 MeV and alter the vacuum sigma mass m a in 
the range 0.6-1.2 GeV because of its uncertainty. The bag 
constant B and dimensionful parameter \o are fixed by 
the vacuum energy density £ = \B = 0.76 GeV fm -3 [15[ 
and the vacuum glueball mass Mq = 1.7 GeV [l6| using 
the following definition: 



Ml 



d 2 V x 
dx 2 



AB 



(2.10) 



The coupling constants Gs and Ga are determined by 
requiring that a nucleon is composed of three constituent 
quarks and a glueball of two constituent gluons, thus, 

M q {T = M = 0) = \m N = 300 MeV, 



M A (T = fi = 0) = -M G 



850 MeV. (2.11) 



3. THERMODYNAMICS 

The model introduced above describes the evolution 
of the two condensates, (a) and (x), driven by temper- 
ature and chemical potential. Figure [T] shows the con- 
tours of the thermodynamic potential, taking the vacuum 
sigma mass being m a — 600 MeV in a-\ plane at ji = 0. 
Increasing temperature from zero, first the system expe- 
riences partial restoration of chiral symmetry at T c hi ra i 
indicated by the dropping a whereas another condensate 
X remains almost a constant. Above T c hirai the potential 
starts to exhibit a meta-stable state at a ~ x ~ and a 
first-order phase transition takes place at T x= o where the 
scale symmetry broken by non-vanishing x is restored. 
Further above this temperature, the system remains at 
the trivial ground state. 

The thermal expectation values of a and x obtained 
from the gap equations in fact show a substantial re- 
duction around the chiral crossover and a jump at the 
first-order transition as seen in Fig. [2] When the sigma 
meson is very massive, A — > oo, one finds 

{<r)*U&, (3.1) 
Xo 

corresponding to non-linear realization of chiral La- 
grangians, and the thermodynamics is governed by a sin- 
gle condensate. Near the chiral symmetry restoration 
point the above relation between the two condensates is 
not expected since the sigma meson cannot be integrated 
out. The condensate of the dilaton field has a weak sen- 
sitivity to a temperature even above the chiral crossover 
and therefore it does not drive the disappearance of the 
chiral condensate. This feature however strongly depends 
on the sigma meson mass and for a larger m CT the gluon 
condensate is more affected by the chiral phase transi- 
tion, as we will discuss below. 

In-medium masses of a and x fields are defined by 

d 2 n . r9 d 2 n 



Mi 



da 2 
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dx 2 
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(3.2) 
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FIG. 1: The contour plots of the thermodynamic potential at finite T and /i = 0: T = 1 53 MeV (chiral crossover), 251 MeV 
(first-order (x) — ► transition) and 300 MeV from left to right. The black circle indicates the ground state. m a — 600 MeV at 
T = was used. 
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FIG. 2: (Color online) The normalized expectation values of 
a and \ 4 fields at fj, = 0. The thick lines are calculated using 
me = 600 MeV at T — and the thin lines using m„ = 900 
MeV. 



Their behavior as functions of temperature is given in 
Fig. [3] Increasing temperature toward T c hirai, M a shows 
a strong sensitivity to the phase transition as observed in 
the standard linear sigma models, whereas M x is rather 
modest. The two masses exhibit a jump when \ vanishes. 
Above this temperature they follow a linear dependence 
of temperature, M a , x ~ T, as expected. 

In Fig. [4] we show the energy density at \i — as 
a function of temperature. The standard linear sigma 
model (LcrM) almost follows the curves below T c hirai, but 
strongly underestimates the Stefan-Boltzmann (SB) limit 
that is a typical drawback of this model. Since the LcrM 
Lagrangian does not contain gluons, bulk thermodynam- 
ics quantities are qualitatively in good agreement with 
the Lattice results when they are normalized by the SB 
limit for massless quarks, whereas not when normalized 
by the SB limit for massless quarks and gluons. What 
we carried out in this paper is to improve the LcrM by in- 
troducing missing gluons. As shown in the figure the SB 
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FIG. 3: (Color online) The thermal masses of a (top) and x 
(bottom) fields at fj, = 0. 



limit is now reproduced. A defect to be removed is too 
strong first-order phase transition even at \i — which is 
absent in Lattice QCD. Also, according to Lattice QCD 
the energy density should approach the SB limit from 
above. We remark that direct comparison must be car- 
ried out in a more realistic framework beyond the mean 
field approximation. As shown in [9(, including thermal 
and quantum fluctuations of meson fields will be partic- 
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FIG. 4: (Color online) The scaled energy density at fj, = 0. 



ularly important around T c . 

The trace anomaly exists at any temperature which is 
the only dimension-full quantity which breaks scale in- 
variance of the theory explicitly. In our model, at high 
temperature and M q <C T the pressure and energy den- 
sity at /i = arc approximately expressed as 



P 
£ 



7tt 2 



1 



-,, T 4 ~ —M%T 2 - -B . 

,q 720 48 9 4 



1 



21^ T 4 _ ll M 2 T 2 1B 

720 48 q 4 



(3.3) 



Consequently, one finds the trace anomaly (interaction 
measure) as 



A(T) 



£ - 3P B i q M t 



(3.4) 



Lattice calculations [13] show that A has a non- 
perturbative term, A ~ 1/T 2 [3l|. We see that this kind 
of contribution comes from the masses of quasiparticles. 
However, the numerical value associated with the effec- 
tive quark mass in Eq. Q3.4p is too small to explain this 
effect. Fluctuations beyond the mean field approxima- 
tion will also contribute to the interaction measure [s| . 

Turning on the quark chemical potential fi practically 
does not affect the temperature at which the gluon con- 
densate vanishes, T x= q 7 whereas the chiral transition 
boundary exhibits an elliptic shape and a critical point 
appears at an intermediate fj,, shown in Fig. [5] (left). The 
boundary line of T x= q in general has a certain /i depen- 
dence via the gap equations. However, the sigma expec- 
tation value above T c hi ra i is small and little affects (x). 
On the other hand, the chiral crossover line gets modi- 
fied significantly depending on m a chosen in vacuum. For 
larger m a the phase boundary is systematically shifted 
to higher T and /x. The critical point also moves toward 
lower T and eventually disappears from the phase dia- 
gram [3]. This is illustrated in Fig. [5] (middle and right). 
The thermodynamics at low temperature and high chem- 
ical potential is essentially same as in the standard linear 
sigma model. 



Making a matching of the trace anomaly between the 
model and QCD would constrain a reliable range of m a . 
The divergence of the dilatation current is given by Q 



= -B 



(X) 
Xo 



T- 



dT 



(3.5) 



The left side of the above equation is mostly saturated 
by the gluon condensate in QCD; 



24 12 ; 



(tW), (3.6) 



where a small contribution due to the explicit breaking 
of chiral symmetry is neglected. Lattice QCD calcula- 
tions show that the thermal gluon condensate decreases 
toward the pseudo-critical temperature of chiral symme- 
try restoration and drops down to a half of its vacuum 
value at T^h-ai, whereas it is quite stable at lower tem- 
peratures [IH . This is also a compatible feature with the 
QCD trace anomaly in terms of the soft and hard dila- 
tons [l9| , i.e. the disappearance of the soft dilaton is as- 
sociated with chiral symmetry restoration and yields the 
melting gluon condensate, or partial restoration of the 
scale symmetry breaking [2fj|. The two equations (|3.5j) 
and (|3.6|) tend to match for a large m a ~ 1 GeV. With 
a small m a the gluon condensate does not show a signif- 
icant drop at T c hirai- Thus, a rather heavy sigma-meson 
in the vacuum seems to be favored by QCD, and this 
is a conceivable scenario known from the vacuum phe- 
nomenology of the scalar mesons. It should be noted that 
the matching is somewhat incomplete; Eq. (|3.5I) exceeds 
Eq. (|3.6[) by ~ 15%. This may indicate that a stronger 
interaction between the quark and gluon sectors should 
be introduced. Besides, updating the gluon condensate 
at finite temperature in Lattice QCD is necessary. 



4. LIMIT OF INFINITELY HEAVY SIGMA 
MESON 

It is instructive to study the phase diagram in the 
A — > oo limit where the sigma meson becomes infinitely 
heavy. As discussed in the previous section, the two 
critical temperatures, T x= o and T^h-ai, get closer with 
increasing m a . With m a ~ 1 GeV they are almost 
on top of each other and larger m a yields an intersec- 
tion of the first-order phase transition of scale symmetry 
and chiral crossover lines at finite /j,. This intersection 
moves to higher /i and lower T for larger m a as shown 
in Fig. [6l The boundary line of scale symmetry restora- 
tion is less sensitive to /i when the chiral symmetry is 
restored. This is because the major /^-dependence comes 
in via the sigma expectation value (a) which is well sup- 
pressed in restored phase. When m a reaches infinity, the 
intersection is kicked out and a single line of the first- 
order phase transition is left. The region where chiral 
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FIG. 5: (Color online) The phase diagram for different vacuum m a : m a = 0.6 GeV (left), 0.9 GeV (middle) and 1.2 GeV 
(right). The filled circle indicates the critical point and the diamond the point where the first-order and crossover lines intersect. 
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FIG. 6: (Color online) The phase diagram for different m^. 
The line notation is same as in Fig. [5] 



symmetry is restored whereas (x) =/= is unfavored in 
this limit. 

The parameters of effective Lagrangians can alter with 
T and \x since they are obtained by integrating higher 
frequency modes out and thus expected to carry infor- 
mation on the underlying QCD. Consequently, the phase 
diagram calculated with the parameters fixed using the 
vacuum quantities would be deformed and the first-order 
phase transition could remain on the phase diagram at 
high /i in a cold system. 



5. IMPLICATIONS FOR THE QCD PHASE 
DIAGRAM 

The present toy model exhibits three regions charac- 
terized by the two condensates: (i) broken phase of chi- 
ral and scale symmetries, (ii) chirally restored but broken 
phase of the scale symmetry because of the non- vanishing 
(x) , and (iii) chirally restored but explicitly broken phase 
of the scale symmetry by temperature. What does the 
thermodynamics of the model suggest concerning the 
QCD phase structure? Vanishing the condensate of the 
dilaton field indicates a disappearance of the gluon com- 



posite at high temperature and its dissociation may sig- 
nal a transition of the system from the confined to decon- 
fined phase. Thus, one identifies the temperature T x= q 
with a temperature at which gluons are released: 

T X=0 ~ T'ifconf • (5.1) 

The model yields a chiral transition temperature that 
is below T x= o i n a wide range of the parameters. In 
Nf = 2 QCD this is compatible with the anomaly match- 
ing which is often used to constrain possible massless ex- 
citations in quantum field theories [21j , and therefore the 
chirally restored phase with confinement is allowed. This 
suggests that the chiral symmetry restoration takes place 
either below or at the deconfinement temperature, i.e. 

Tdiiral ~ ^deconf ' 0>>2) 

where at T^ coai the quarks are released whereas the glu- 
ons remain confined and it is not necessarily equal to 
-^deconf ^ s we nave seen in the previous section, a large 
m a can match with the QCD requirement at \i = 0. 
This leads to the three distinct temperatures which may 
be close to each other on the phase diagram. We note 
that this is consistent with the recent observation using 
a renormalization group analysis where the fixed point 
of four-fermion interactions associated with confinement 
plays an essential role [22j. At finite /i no reliable con- 
straint from QCD is known. A suggestive phase diagram 
is given in Fig. [7J 

6. CONCLUSIONS AND REMARKS 

In this paper we have studied thermodynamics and 
the phase structure of a QCD-like model whose degrees 
of freedom are constituent quarks and gluons. Both chi- 
ral and scale symmetries are implemented in the model 
by introducing mean fields representing qq and G^ U G^ V . 
These symmetries are dynamically broken at low temper- 
ature and density. The model thus mimics the features 
of QCD in the strong coupling region, i.e. the sponta- 
neous breaking of chiral symmetry and trace anomaly. 
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FIG. 7: (Color online) A sketch of the QCD phase diagram. 



The results suggest that a system in deconfined phase 
develops gradually with increasing temperature/density 
toward weakly-interacting quark-gluon matter composed 
of almost massless quarks and gluons. 

The condensates of the sigma and dilaton fields are 
dynamically linked via their gap equations. How strong 
they are correlated depends crucially on the sigma-meson 
mass m a chosen in vacuum. We found that a large 
m a ~ 1 GeV is consistent with the lattice result re- 
garding the thermal behavior of the gluon condensate. 
This further leads to the chiral phase transition which 
takes place almost simultaneously with the deconfinc- 
ment transition at fi ~ 0. At finite fi these two tran- 
sitions are expected to be separated. 

In the scalar sector of low-mass hadrons, scalar quarko- 
nium, tetra-quark states [23[ and glueballs are expected 
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